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We combine results from CDF and DO's direct searches for the standard model (SM) Higgs boson 
{H) produced in pp collisions at the Fermilab Tevatron at y/s = 1.96 TeV. The results presented 
here include those channels which are most sensitive to Higgs bosons with mass between 130 and 
200 GeV/c^, namely searches targeted at Higgs boson decays to W^W~ , although acceptance for 
decays into r^r~ and 77 is included. Compared to the previous Tevatron Higgs search combination, 
more data have been added and the analyses have been improved to gain sensitivity. We use the 
MSTW08 parton distribution functions and the latest gg ^ H theoretical cross section predictions 
when testing for the presence of a SM Higgs boson. With up to 7.1 fb~^ of data analyzed at CDF, 
and up to 8.2 fb^^ at DO, the 95% C.L. upper limits on Higgs boson production is a factor of 0.54 
times the SM cross section for a Higgs boson mass of 165 GeV/c^ We exclude at the 95% C.L. the 
region 158 < mn < 173 GeV/c^ 



Preliminary Results 



* The Tevatron New-Phenomena and Higgs Working Group can be contacted at TEVNPHWG@fnal.gov. More information can be found 
at http://tevnphwg.fnal.gov/ 
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I. INTRODUCTION 

The search for a mechanism for electroweak symmetry breaking, and in particular for a standard model (SM) Higgs 
boson, has been a major goal of particle physics for many years, and is a central part of the Fermilab Tevatron physics 
program, and also that of the Large Hadron Collider. Recently, the ATLAS and CMS collaborations at the Large 
Hadron Collider have released results on searches for the Standard Model Higgs boson decaying to W^W" [5), Ig], 
ZZ J] , 77 [8] Both the CDF and DO collaborations have performed new combinations [H l2l| of multiple direct searches 
for the SM Higgs boson. The new searches cover a larger data sample, and incorporate improved analysis techniques 
compared to previous analyses. The sensitivities of these new combinations significantly exceed those of previous 
combinations (3), |^] . 

In this note, we combine the most recent results of searches for the SM Higgs boson produced in pp collisions 
at y/s ~ 1.96 TeV. The analyses combined here seek signals of Higgs bosons produced in association with vector 
bosons {qq — > W/ZH), through gluon-gluon fusion {gg — ^ H), and through vector boson fusion (VBF) {qq — >■ q'q'H) 
corresponding to integrated luminosities up to 7.1 fb~^ at CDF and up to 8.2 fb^^ at DO. The searches included 
here target Higgs boson decays to W~^W~ , although acceptance for decays into r+r^ and 77 are included in the DO 
channels. 

The searches are separated into 46 mutually exclusive final states (12 for CDF and 34 for DO; see Tables III and IhI) 
referred to as "analysis sub-channels" in this note. The selection procedures for each analysis are detailed in Refs. ji| 
through [3], and are briefly described below. 

II. CHANNELS INCLUDED IN THE COMBINATION 

Event selections are similar for the corresponding CDF and DO analyses. 

For the H — > W^W^ analyses which seek events in which both W bosons decay leptonically, signal events are 
characterized by large l^rp and two oppositely-signed, isolated leptons. The presence of neutrinos in the final state 
prevents the accurate reconstruction of the candidate Higgs boson mass. DO selects events containing large and 
electrons and/or muons, dividing the data sample into three final states: e+e~, e^/i^, and fj,~^fi~. The searches in 
the three final states each use 8.1 fb~^ of data which is subdivided according to the number of jets in the event: 0, 
1, or 2 or more ("2-I-") jets. Decays involving tau leptons are included in two orthogonal ways. For the first time, 
a dedicated analysis {^Thad) using 7.3 fb^^ of data studying the final state involving a muon and a hadronic tau 
decay, is included in the Tevatron combination. Final states involving other tau decays and mis-identified hadronic 
tau decays are included in the e+e^, e^/i^, and fi^id^ final state analyses. CDF separates the H — > W^W^ events in 
five non-overlapping samples, split into both "high s/fo" and "low s/&" categories based on lepton types and different 
categories based on the number of reconstructed jets: 0, 1, or 2+ jets. The sample with two or more jets is not 
split into low s/b and high s/b lepton categories due to low statistics. A sixth CDF channel is the low dilepton 
mass (m^+i-) channel, which accepts events with mg+g- < 16 GeV. This channel increases the sensitivity of the 
H — >• W^W~ analyses at low m//, adding 10% additional acceptance at mn — 120 GeV. The division of events into 
jet categories allows the analysis discriminants to separate three different categories of signals from the backgrounds 
more effectively. The signal production mechanisms considered are gg ^ H W^W~, WH + ZH — ^ jjW^W", 
and vector-boson fusion. The DO e+e~, e^fjr^, and /x^/x" final state channels use boosted decision tree outputs as 
the final discriminants while the pLThad channel uses neural networks. CDF uses neural-network outputs, including 
likelihoods constructed from calculated matrix-element probabilities as additional inputs for the 0-jet bin. 

DO includes VH — -I- AT analyses in which the associated vector boson and the W boson from the Higgs boson 
decay that has the same charge are required to decay leptonically, thereby defining three like-sign dilepton final states 
(e^e^, e^/x^, and /x^/i^). The combined output of two decision trees, trained against the instrumental and diboson 
backgrounds respectively, is used as the final discriminant. CDF also includes a separate analysis of events with 
same-sign leptons and large to incorporate additional potential signal from associated production events in which 
the two leptons (one from the associated vector boson and one from a W boson produced in the Higgs decay) have 
the same charge. CDF for the first time also incorporates three tri-lepton channels to include additional associated 
production contributions where leptons result from the associated W boson and the two W bosons produced in the 
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Higgs decay or where an associated Z boson decays into a dilepton pair and a third lepton is produced in the decay 
of either of the W bosons resulting from the Higgs decay. In the latter case, CDF separates the sample into 1 jet and 
2+ jet sub-channels to fully take advantage of the Higgs mass constraint available in the 2+ jet case where all of the 
decay products are reconstructed. 

CDF also includes opposite-sign channels in which one of the two lepton candidates is a hadronic tau. Events 
are separated into e-r and [i-t channels. The final discriminants are obtained from boosted decision trees which 
incorporate both hadronic tau identification and kinematic event variables as inputs. 

DO includes channels in which one of the W bosons in the R — > W^W^ process decays leptonically and the other 
decays hadronically. Electron and muon final states are studied separately, each with 5.4 fb"-'^ of data. Random 
forests are used for the final discriminants. 

DO has updated its search for Higgs boson production in which the Higgs decays into a pair of photons to include 
8.2 fb~^ of data and to use a boosted decision tree as the final discriminant. 

The DO analysis of the r+r" -I- 2 jets final state has also been updated to include additional data (4.3 fb^^) and 
both muonic and electronic tau decays. The output of boosted decision trees is used as the final discriminant. 

For both CDF and DO, events from QCD multijet (instrumental) backgrounds are measured in independent data 
samples using different methods. For CDF, backgrounds from SM processes with electroweak gauge bosons or top 
quarks were generated using pythia ALPGEN [H], mc@nlo \Vh. and herwig [i1| programs. For DO, these 



backgrounds were generated using pythia, ALPGEN, and COMPHEP [19|, with pythia providing parton-showering 
and hadronization for all the generators. These background processes were normalized using either experimental data 
or next-to- leading order calculations (including mcfm f23| for the W+ heavy flavor process). 

Tables [11 and HIl summarize, for CDF and DO respectively, the integrated luminosities, the Higgs boson mass ranges 
over which the searches are performed, and references to further details for each analysis. 



TABLE I: Luminosity, explored mass range and references for the different processes and final states {J. = e, ^i) for the CDF 
analyses. 



Channel 


Luminosity (fb ^) 


tuh range (GeV/c^) 


Reference 


H W^W 2x (0,1 jets)-h(2-h jets) + (low-m^^) + (e-r,w) + (M-T-had) 


7.1 


130-200 


yj 


WH — ^ WW~^W~ (same-sign leptons 1+ jets)-l-(tri-leptons) 


7.1 


130-200 




ZH ZW+W~ (tri-leptons 1 jet) + (tri-leptons 2+ jets) 


7.1 


130-200 





TABLE IL Luminosity, explored mass range and references for the difi'erent processes and final states {£ = e, fi) for the DO 
analyses. Most analyses are in addition analyzed separately for Run Ila and lib. 



Channel Luminosity (fb ^) iriH range (GeV/c^) Reference 

H -> W+W- t^uf^v (0,1,2-F jet) O 130-200 

H W+W- -f ^LVTHadV 7.3 130-200 [10] 

H ^W+W- ^ Ivjj 5.4 130-200 [llj 

VH-fl+l+ +X 5.3 130-200 [12] 

H+X^e+r^^Jj 4.3 130-200 [13] 

i/ ^ 77 8.2 130-150 [IJ 
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III. SIGNAL PREDICTIONS AND UNCERTAINTIES 



We normalize our Higgs boson signal predictions to the most recent highest-order calculations available, for all 
production processes considered. The largest production cross section, a^gg — >■ H), is calculated at next-to-next-to- 
leading order (NNLO) in QCD with soft gluon resummation to next-to-next-to-leading-log (NNLL) accuracy, and 
also includes two- loop electroweak effects and handling of the running b quark mass [2l|, ■ The numerical values 
in Table Hill are updates [1^ of these predictions with rrit set to 173.1 GeV/c^ [13], and an exact treatment of the 
massive top and bottom loop corrections up to next-to-leading-order (NLO) + next-to-leading-log (NLL) accuracy. 
The factorization and renormalization scale choice for this calculation is fj,p = IJ'R = rnu- These calculations are 
refinements of earlier NNLO calculations of the gg ^ H production cross section [25l - l27| . Electroweak corrections 
were computed in Refs. [H, 111]. Soft gluon resummation was introduced in the prediction of the gg H production 
cross section in Ref. "s^. 

The gg ^ H production cross section depends strongly on the gluon parton density function, and the accompanying 
value of as (9^ ) • The cross sections used here are calculated with the MSTW 2008 NNLO PDF set [Slj , as recommended 
by the PDF4LHC working group f32l|. The reason to use this PDF parameterization is that it is the only NNLO 
prescription that results from a fully global fit to all relevant data. We follow the PDF4LHC working group's 
prescription to evaluate the uncertainties on the gg H production cross section due to the PDFs. This prescription 
is to evaluate the predictions of a{gg H) at NLO using the global NLO PDF sets CTEQ6.6 ^33)], MSTW08 [sll, 
and NNPDF2.0 [3J|, and to take the envelope of the predictions and their uncertainties due to PDF-|-as for the 
three sets as the uncertainty range at NLO. The ratio of the NLO uncertainty range to that of just MSTW08 is then 
used to scale the NNLO MSTW08 PDF-l-as uncertainty, to estimate a larger uncertainty at NNLO. This procedure 
roughly doubles the PDF-has uncertainty from MSTW08 at NNLO alone. Other PDF sets, notably ABKM09 [Hi 
and HERAPDFl.O ^i^, have been proposed [s^, as alternate PDF sets for use in setting the uncertainties on 
(j{gg ^ H). These two PDF sets fit HERA data only in the case of HERAPDFl.O, and HERA, fixed-target deep- 
inelastic scattering (DIS), fixed-target Drell-Yan (DY), and neutrino- nucleon di-muon data in the case of ABKM09, 
and do not include Tevatron jet data [s^jlll. We note that HERA, fixed-target DIS, fixed-target DY, and neutrino- 
nucleon di-muon data do not probe directly the gluon PDF at the scales most relevant for our analyses, but are 
most sensitive to the quark distributions. The high-i?T Tevatron jet data contain gluon-initiated processes at the 
relevant momentum transfer scales. The PDF constraints from the HERA, fixed-target DIS, fixed-target DY, and 
di-muon data are extrapolated to make predictions of high-x gluons at the relevant energy scales. The reliability of 
these extrapolations is checked for the HERAPDFl.O PDF set by showing the NLO prediction 42] of the inclusive jet 
spectrum as a function of jet Et and ry, as compared with measurements by CDF and DO as shown in Figures [T] and [H 
respectively. The measured inclusive jet cross sections are significantly under- predicted using this PDF set. The 
corresponding comparisons using the same NLO jet cross section calculation with the MSTW08 PDF set are shown 



in Figures 49 and 50 in Ref. |3ll- The agreement with the jet data is much better with the MSTW2008 PDF set. 



which is expected, since those data are included in the MSTW08 PDF fit. Consequently the high- a; gluon component 
is larger in the MSTW2008 PDF set than in the HERAPDFl.O PDF set. This forms an important experimental 
foundation for the prediction of the giuon-fusion Higgs boson production cross section. Recently, it has been pointed 
out that the gluon PDF in the ABKM09 fit, and thus the values of (j{gg — i/), depend on the handling of the 
data input from the New Muon Collaboration (NMC), specifically whether the differential spectrum as measured by 
NMC is included in the fit or whether the structure function F2 reported by NMC is used in the fit. It is verified by 
the NNPDF Collaboration Q that the NNPDF2.1 global PDF fit is negligibly impacted by changing the handling of 
the input NMC data, indicating that the global fit, which includes Tevatron jet data, is more solid and less sensitive 
to the precise handling of NMC data. It should be noted however that this study is performed only at NLO. 

We also include uncertainties on u{gg — >■ H) due to uncalculated higher order processes by following the standard 
procedure of varying the factorization renormalization scales up and down together by a factor k = 2. It had been 
proposed to vary /i/j and independently (STj . but the maximum impact on (T{gg — s- H) is found when varying the 
two scales together. At a central scale choice of = = mH/2, the fixed-order calculations at NLO and NNLO 
converge rapidly and the scale uncertainties generously cover the differences. Furthermore, the NNLO calculation 
using a scale choice of fin — fip — m^f/2 is very close to the NNLL-I-NNLO calculation that we use. The authors of 
the NNLL-I-NNLO calculation recommend that we use the scale choice = ^p = niH for the central value [Hj. The 
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uncertainty on <7{gg — >■ H) is evaluated by varying f^n — from a lower value of mull to an upper value of IrriH^ 
the customary variation, and by evaluating the NNLO+NNLL cross sections at these scales. A larger scale variation 
has been proposed, varying the scale by a factor k = 3 [s^l- The original justification for this proposal [s^l was for the 
LO calculation to cover the NNLO central value. The LO calculation however does not have enough scale-dependent 
terms for its scale uncertainty to be fairly considered - it is missing significant QCD radiation and other terms that 
increase the scale dependence at NLO compared with that at LO. The second justification [ssj is simply that the 
higher-order QCD corrections are larger than those of other processes, such as Drell-Yan production. The corrections 
are large due to the fact that there are gluons in the initial state which radiate more gluons more frequently than 
the quarks that initiate other processes. The scale dependence of the calculations of <j{gg — >■ H) which include this 
radiation and other terms is correspondingly larger, and the customary range of scale variation provides uncertainties 
which are commensurate with the values of the corrections. An independent calculation by Ahrens, Becher, Neubert, 
and Yang [4^ using a renormalization-group improved resummation of radiative corrections at NNLL accuracy gives 
an alternative calculation of <j{gg — !■ i?) that converges more rapidly and requires even smaller scale uncertainties. 
The central values of this calculation of ai^gg — !■ i?) are in excellent agreement with those used in our analysis and 
which are given in Table IIIII 

Because the H — > W'^W^ analyses separate the data into categories based on the numbers of observed jets, we 
assess factorization and renormalization scale and PDF-l-as variations separately for each jet category as evaluated 
in Ref. [46|]. This calculation is at NNLO for H+0 jets, at NLO for H+1 jet, and at LO for H+2 or more jets. A 
newer, more precise calculation [43| of the H+2 or more jets cross section at NLO is used to evaluate the uncertainties 
in this category. These scale uncertainties are used instead of the inclusive NNLL+NNLO scale uncertainty because 
we require them in each jet category, and the uncertainties we use are significantly larger than the inclusive scale 
uncertainty. The scale choice affects the pr spectrum of the Higgs boson when produced in gluon-gluon fusion, 
and this effect changes the acceptance of the selection requirements and also the shapes of the distributions of the 
final discriminants. The effect of the acceptance change is included in the calculations of Ref. [1^ and Ref. [13], as 
the experimental requirements are simulated in these calculations. The effects on the final discriminant shapes are 
obtained by reweighting the px spectrum of the Higgs boson production in our Monte Carlo simulation to higher-order 



calculations. The Monte Carlo signal simulation used by CDF and DO is provided by the LO generator pythia [15 | 
which includes a parton shower and fragmentation and hadronization models. We reweight the Higgs boson pT 
spectra in our pythia Monte Carlo samples to that predicted by hqt 48] when making predictions of differential 
distributions of ag — » H signal events. To evaluate the impact of the scale uncertainty on our differential spectra, we 
use the resbos [49] generator, and apply the scale-dependent differences in the Higgs boson px spectrum to the hqt 
prediction, and propagate these to our final discriminants as a systematic uncertainty on the shape, which is included 
in the calculation of the limits. 

We treat the scale uncertainties as 100% correlated between jet categories and between CDF and DO, and also treat 
the PDF-|-as uncertainties in the cross section as correlated between jet categories and between CDF and DO. We 
treat however the PDF-|-as uncertainty as uncorrelated with the scale uncertainty. The main reason is that the PDF 
uncertainty arises from experimental uncertainties and PDF parameterization choices, while the scale uncertainty 
arises from neglected higher-order terms in the perturbative cross section calculations. The PDF predictions do 
depend on the scale choice, however. The dependence of the prediction of (T{gg — ^ H) due to the scale via the 
PDF dependence is included as part of the scale uncertainty and not as part of the PDF uncertainty, to ensure that 
all scale dependence is considered correlated. Furthermore, we have verified [sO] that the fractional change in the 
prediction of (j{gg — H) due to the PDF variation depends negligibly on the value of the scale choice, justifying 
the treatment of PDF and scale uncertainties as uncorrelated. As described in Section |Vl systematic uncertainties 
arising from uncorrelated sources should not be added linearly together, but instead should be considered to fluctuate 
independently of one another. We have investigated the use of a Gaussian prior for the scale uncertainty using the 
K = 2 variation as the ±la variation, and compared it with the use of a flat prior between the low edge and the upper 
edge, and have found a negligible impact on our sensitivity and our observed limits. The main reason for this is that 
while the flat prior has a higher density near the edges, the Gaussian prior has longer tails, and both effects essentially 
compensate. The use of a flat prior with edges for the scale uncertainty does not affect the choice of whether to add 
this uncertainty linearly or not with the PDF uncertainty; that choice is determined by the correlations of the two 
sources of uncertainty. 

Another source of uncertainty in the prediction of (j{gg — )■ H) is the extrapolation of the QCD corrections computed 
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TABLE III: The production cross sections and decay branching fractions for the SM Higgs boson assumed for the combination. 



rriH 




cywH 


CrZH 


0"VSF 


B{H T+T-) 


B(H W+W-) 


B{H ZZ) 


B{H -5> 77) 


(GeV/c^) 


(fb) 


(fb) 


(fb) 


(fb) 


(%) 


(%) 


(%) 


(%) 


130 


842.9 


112.00 


68.5 


62.1 


5.305 


29.43 


3.858 


0.2182 


135 


750.8 


97.20 


60.0 


57.5 


4.400 


39.10 


5.319 


0.2077 


140 


670.6 


84.60 


52.7 


53.2 


3.472 


49.16 


6.715 


0.1897 


145 


600.6 


73.70 


46.3 


49.4 


2.585 


59.15 


7.771 


0.1653 


150 


539.1 


64.40 


40.8 


45.8 


1.778 


68.91 


8.143 


0.1357 


155 


484.0 


56.20 


35.9 


42.4 


1.057 


78.92 


7.297 


0.09997 


160 


432.3 


48.50 


31.4 


39.4 


0.403 


90.48 


4.185 


0.05365 


165 


383.7 


43.60 


28.4 


36.6 


0.140 


95.91 


2.216 


0.02330 


170 


344.0 


38.50 


25.3 


34.0 


0.093 


96.39 


2.351 


0.01598 


175 


309.7 


34.00 


22.5 


31.6 


0.073 


95.81 


3.204 


0.01236 


180 


279.2 


30.10 


20.0 


29.4 


0.059 


93.25 


5.937 


0.01024 


185 


252.1 


26.90 


17.9 


27.3 


0.046 


84.50 


14.86 


0.008128 


190 


228.0 


24.00 


16.1 


25.4 


0.038 


78.70 


20.77 


0.006774 


195 


207.2 


21.40 


14.4 


23.7 


0.033 


75.88 


23.66 


0.005919 


200 


189.1 


19.10 


13.0 


22.0 


0.029 


74.26 


25.33 


0.005285 



for the heavy top quark loops to the hght-quark loops included as part of the electroweak corrections. Uncertainties 



at the level of 1-2% are already included in the cross section values we use [2l|, [2^. In Ref. [2l|, it is argued that the 
factorization of QCD corrections is known to work well for Higgs boson masses many times in excess of the masses of 
the loop particles. A 4% change in the predicted cross section is seen when all QCD corrections are removed from the 
diagrams containing light-flavored quark loops, which is too conservative. For the h quark loop, which is computed 
separately in Ref. the QCD corrections are much smaller than for the top loop, further giving confidence that it 
does not introduce large uncertainties. 

We include all significant Higgs boson production modes in our searches. Besides gluon-gluon fusion through 
virtual quark loops (GGF), we include Higgs boson production in association with a W or Z vector boson, and 
vector boson fusion (VBF). We use the WH and ZH production cross sections computed at NNLO in Ref. [sH . 
This calculation starts with the NLO calculation of v2hv [s^ and includes NNLO QCD contributions [s^, as well 
as one- loop electroweak corrections [H^. We use the VBF cross section computed at NNLO in QCD in Ref. jssj . 
Electroweak corrections to the VBF production cross section are computed with the hawk program [s^, and are 
small and negative (-2% to -3%) for the Higgs boson mass range considered here. They are not included in this 
update but will be incorporated in future results. 

In order to predict the distributions of the kinematics of Higgs boson signal events, CDF and DO use the pythia fllj 
Monte Carlo program, with CTEQ5L and CTEQ6L ^] leading-order (LO) parton distribution functions. The Higgs 
boson decay branching ratio predictions are calculated with HDECAY [53| , and are also listed in Table IIIIl We use 
HDECAY Version 3.53. While the HWW coupling is well predicted, Br{H — >■ depends on the partial widths 

of all other Higgs boson decays. The partial width r{H — 66) is sensitive to rrib and as, T{H — >■ cc) is sensitive 
to rric and as, and T{H — )■ gg) is sensitive to ag. The impacts of these uncertainties on Br{H — W~^W~) depend 
on rriH due to the fact that Br{H — )■ 66), Br{H — cc), Br{H — )■ gg) become very small for Higgs boson masses 
above 160 GeV/c^, while they have a larger impact for lower mn- We use the uncertainties on the branching fraction 
Br{H — ^ VF+VF~) from Ref. At niH = 130 GeV/c^, for example, the mt variation gives a ^i;7o% relative 

variation in Br{H — >• W'^W~), as gives a +i;o9% variation, and rric gives a +o;5i% variation. At mn — 165 GeV/c^, 
all three of these uncertainties are below 0.1% and remain small for higher tuh- 
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Tevatron Jet Cross Sections 
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P^J^^' [GeV/c] 



FIG. 1: Inclusive jet cross section measurements by DO as a function of jet pt in six bins of jet rapidity y, divided by the 
theoretical prediction. The data systematic uncertainties are displayed by the full shaded band. The theoretical prediction is 
an NLO perturbative QCD calculation with renormalization and factorization scales set to jet pr using the HERAPDFl.O 
PDFs [3a][40l] and including nonperturbative corrections. The PDF uncertainties are shown as hatched regions - the blue 
hatching shows the total PDF uncertainty, and the red hatching shows the portion of the total due to experimental uncertainties. 
This figure is obtained using the same data and methods described in Ref. [iJl . 
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Tevatron Jet Cross Sections 
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FIG. 2: Inclusive jet cross section measurements by CDF as a function of jet pr in six bins of jet rapidity y, divided by 
the theoretical prediction. The theoretical prediction is an NLO perturbative QCD calculation [i^], with renormalization and 
factorization scales set to jet pr using the HERAPDFl.O PDFs [3^[43] and including nonperturbative corrections. The PDF 
uncertainties are shown as hatched regions - the blue hatching shows the total PDF uncertainty, and the red hatching shows 
the portion of the total due to experimental uncertainties. The CDF data are published in Ref. 39]. Experimental uncertainties 
on the CDF data are comparable to those of the DO data, shown in Fig. [1] 
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Tevatron Run II Preliminary, L < 8.2 tb ' 
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FIG. 3: Distributions of logjQ(s/6), for the data from all contributing channels from CDF and DO, for Higgs boson masses of 
150 and 165 GeV/c^ The data are shown with points, and the expected signal is shown stacked on top of the backgrounds. 
Underflows and overflows are collected into the bottom and top bins. 



IV. DISTRIBUTIONS OF CANDIDATES 



All analyses provide binned histograms of the final discriminant variables for the signal and background predictions, 
itemized separately for each source, and the observed data. The number of channels combined is large, and the number 
of bins in each channel is large. Therefore, the task of assembling histograms and checking whether the expected and 
observed limits are consistent with the input predictions and observed data is difficult. We thus provide histograms 
that aggregate all channels' signal, background, and data together. In order to preserve most of the sensitivity gain 
that is achieved by binning the data instead of collecting them all together and counting, we aggregate the data and 
predictions in narrow bins of signal-to-background ratio, s/b. Data with similar s/b may be added together with 
no loss in sensitivity, assuming similar systematic errors on the predictions. The aggregate histograms do not show 
the effects of systematic uncertainties, but instead compare the data with the central predictions supplied by each 
analysis. 

The range of s/b is quite large in each analysis, and so log]^g(s/6) is chosen as the plotting variable. Plots of the 
distributions of logio(s/&) are shown for rriH = 150 and 165 GeV/c^ in Figure[3l These distributions can be integrated 
from the high-s/6 side downwards, showing the sums of signal, background, and data for the most pure portions of the 
selection of all channels added together. These integrals can be seen in Figure |H The most significant candidates are 
found in the bins with the highest s/b; an excess in these bins relative to the background prediction drives the Higgs 
boson cross section limit upwards, while a deficit drives it downwards. The lower-s/5 bins show that the modeling of 
the rates and kinematic distributions of the backgrounds is very good. 

We also show the distributions of the data after subtracting the expected background, and compare that with the 
expected signal yield for a standard model Higgs boson, after collecting all bins in all channels sorted by s/b. These 
background-subtracted distributions are shown in Figure [51 These graphs also show the remaining uncertainty on the 
background prediction after fitting the background model to the data within the systematic uncertainties on the rates 
and shapes in each contributing channel's templates. 
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FIG. 4: Integrated distributions of s/b, starting at the liigii s/h side, for Higgs boson masses of 150 and 165 GeV/c^ The total 
signal-plus-background and background-only integrals are shown separately, along with the data sums. Data are only shown 
for bins that have data events in them. 



1 11 Preliminary-4- Data-Background' 
I I SM Higgs Signal 
±1 s.d. on Background 




0.6 -0.4 -0.2 
logjj(s/b) 



m 


50 


o 


40 




c 


30 


> 


20 




10 









-10 




-20 




-30 




-40 




-50 



Tevatron Run II Preliminar 
L < 8.2 tV' 



f-\— Data-Background 
I I SM Higgs Signal 

±1 s.d. on Background 




March 7,2011 



mH=165 GeWlc 



-1.6 -1.4 -1.2 -1 -O.i 



-0.6 -0.4 -0.2 
logio(S' 



FIG. 5: Background-subtracted data distributions for all channels, summed in bins of s/b, for Higgs boson masses of 150, 
and 165 GeY/c?. The background has been fit, within its systematic uncertainties, to the data. The points with error bars 
indicate the background-subtracted data; the sizes of the error bars are the square roots of the predicted background in each 
bin. The unshaded (blue-outline) histogram shows the systematic uncertainty on the best-fit background model, and the shaded 
histogram shows the expected signal for a Standard Model Higgs boson. 



V. COMBINING CHANNELS 



To gain confidence that the final result does not depend on the details of the statistical formulation, we perform 
two types of combinations, using Bayesian and Modified Frequentist approaches, which yield limits on the Higgs 
boson production rate that agree within 10% at each value of ran, and within 1% on average. Both methods rely on 
distributions in the final discriminants, and not just on their single integrated values. Systematic uncertainties enter 
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on the predicted number of signal and background events as well as on the distribution of the discriminants in each 
analysis ("shape uncertainties"). Both methods use likelihood calculations based on Poisson probabilities. 

Both methods treat the systematic uncertainties in a Bayesian fashion, assigning a prior distribution to each 
source of uncertainty, parameterized by a nuisance parameter, and propagating the impacts of varying each nuisance 
parameter to the signal and background predictions, with all correlations included. A single nuisance parameter 
may affect the signal and background predictions in many bins of many different analyses. Independent nuisance 
parameters are allowed to vary separately within their prior distributions. Both methods use the data to constrain 
the values of the nuisance parameters, one by integration, the other by fitting. These methods reduce the impact 
of prior uncertainty in the nuisance parameters thus improving the sensitivity. Because of these constraints to the 
data, it is important to evaluate the uncertainties and correlations properly, and to allow independent parameters to 
vary separately, otherwise a fit may overconstrain a parameter and extrapolate its use improperly. The impacts of 
correlated uncertainties add together linearly on a particular prediction, while those of uncorrelated uncertainties are 
convoluted together, which is similar to adding in quadrature. Adding uncorrelated uncertainties linearly implies a 
correlation which is not present and which may result in incorrect results. 

Both methods set limits at the 95% confidence (or "credibility" in the case of the Bayesian method). This is 
a probabalistic statement. To be consistent and accurate, we give equal treatment to all sources of uncertainty, 
both theoretical and experimental. Ref. (ssj proposes testing the minimum possible production cross section and 
quoting our exclusion limits based on only this prediction. We do not claim that our exclusion limits are independent 
of the model choices that are made, and hence we are not required to only test the case in which all values of 
all uncertain parameters fluctuate simultaneously to the values that produce the weakest results. To do so would 
provide inconsistent results when computing limits and when attempting to discover. A set of values for the nuisance 
parameters which weakens the limit may strengthen a discovery and vice versa, particularly those parameters that 
affect the background rate and shape predictions. By always setting uncertain parameters to their most extreme 
values, we could find we have an excess of data over the background prediction when we set our limits, and a deficit 
of data with respect to the background in the same sample when attempting to make a discovery. The prescriptions 
described below provide a consistent method for both tasks. 



A. Bayesian Method 

Because there is no experimental information on the production cross section for the Higgs boson, in the Bayesian 
technique QjlOl we assign a flat prior for the total number of selected Higgs events. For a given Higgs boson mass, the 
combined likelihood is a product of likelihoods for the individual channels, each of which is a product over histogram 
bins: 

Nc Ni "np 

C{R, s, b\n, e) X -,{9) = n n X n (1) 

1=1 j=l fe=l 

where the first product is over the number of channels [Nc), and the second product is over Ni, histogram bins 
containing events, binned in ranges of the final discriminants used for individual analyses, such as the dijet mass, 
neural-network outputs, or matrix-element likelihoods. The parameters that contribute to the expected bin contents 
are ^ij = R x Sij{0) -\- bij{9) for the channel i and the histogram bin j, where Sij and hij represent the expected 
background and signal in the bin, and i? is a scaling factor applied to the signal to test the sensitivity level of the 
experiment. Truncated Gaussian priors are used for each of the nuisance parameters Ok, which define the sensitivity of 
the predicted signal and background estimates to systematic uncertainties. These can take the form of uncertainties 
on overall rates, as well as the shapes of the distributions used for combination. These systematic uncertainties can 
be far larger than the expected SM Higgs boson signal, and are therefore important in the calculation of limits. The 
truncation is applied so that no prediction of any signal or background in any bin is negative. The posterior density 
function is then integrated over all parameters (including correlations) except for R, and a 95% credibility level upper 
limit on R is estimated by calculating the value of R that corresponds to 95% of the area of the resulting distribution. 
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B. Modified Frequentist Metliod 



The Modified Frequentist technique relies on the CLg method, using a log-likelihood ratio (LLR) as test statistic [2|: 




where Hi denotes the test hypothesis, which admits the presence of SM backgrounds and a Higgs boson signal, while 
Ho is the null hypothesis, for only SM backgrounds. The probabilities p are computed using the best-fit values of the 
nuisance parameters for each pseudo-experiment, separately for each of the two hypotheses, and include the Poisson 
probabilities of observing the data multiplied by Gaussian priors for the values of the nuisance parameters. This 
technique extends the LEP procedure |60j which does not involve a fit, in order to yield better sensitivity when 
expected signals are small and systematic uncertainties on backgrounds are large [Blf. 

The CLg technique involves computing two p- values, CLg+b and CLb. The latter is defined by 



where LLi?obs is the value of the test statistic computed for the data. 1 — CLb is the probability of observing a 
signal-plus-background-like outcome without the presence of signal, i.e. the probability that an upward fluctuation of 
the background provides a signal-plus-background-like response as observed in data. The other p-value is defined by 



and this corresponds to the probability of a downward fluctuation of the sum of signal and background in the data. 
A small value of CLg+b reflects inconsistency with Hi. It is also possible to have a downward fiuctuation in data 
even in the absence of any signal, and a small value of CLg+b is possible even if the expected signal is so small 
that it cannot be tested with the experiment. To minimize the possibility of excluding a signal to which there is 
insufficient sensitivity (an outcome expected 5% of the time at the 95% C.L., for full coverage), we use the quantity 
CLs = CLs+b/CLb. If CLs < 0.05 for a particular choice of Hi, that hypothesis is deemed to be excluded at the 
95% C.L. In an analogous way, the expected CLb, CLg+b and CLg values are computed from the median of the LLR 
distribution for the background-only hypothesis. 

Dividing CLg+b by CLb incurs a median penalty of a factor of two in the expected p- value, as the distribution of CLb 
is uniform under the null hypothesis. One may use CLg+b directly to set limits, and solve the problem of excluding 
non-testable models by introducing a "power-constraint". The power-constrained limit (PCL) method [1,0,[1] involves 
excluding models for which CLg+b < 0.05, but if the limit thus obtained makes an excursion below a previously chosen 
lower bound, to quote the lower-bound limit instead. The ATLAS collaboration sets this lower bound at —la in the 
distribution of limits expected from background-only outcomes. The PCL method thus retains the desired coverage 
rate and does not exclude untestable models, and it provides by construction stronger expected and observed limits 
than CLg. The expected and observed PCL limits with Tevatron data may be estimated from FigurelHland Table IVllIl 
which show and list the observed and expected values of CLs_|_b as functions of mn- Nowhere in the range is there 
more than a la downward fluctuation relative to the background prediction, and so a similar power constraint to 
the one ATLAS applies would not have an effect on the observed limits. With the PCL method, the region of mn 
expected to be excluded at the 95% C.L. grows by ^40% compared to that obtained with the CLg method. The PCL 
method has a larger false exclusion rate than the CLg and Bayesian methods that we use to quote our results, the 
distribution of possible limits under the null hypothesis is highly asymmetrical, and this distribution depends on the 
arbitrary choice of the location of the power constraint. We decided to continue quoting CLg and Bayesian limits also 
because they provide a strong numerical cross-check of each other, and they are comparable to those of LEP [6^ . 
CMS @, and our own previous results 01 ■ 

Systematic uncertainties are included by fluctuating the predictions for signal and background rates in each bin of 
each histogram in a correlated way when generating the pseudo-experiments used to compute CLg+b and CLb. 



1 - CLb = p{LLR < LLR,^,\Ho), 



(3) 



CLs+b ^ p{LLR> LLR,^,\Hi), 



(4) 
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C. Systematic Uncertainties 

Systematic uncertainties differ between experiments and analyses, and they affect the rates and shapes of the 
predicted signal and background in correlated ways. The combined results incorporate the sensitivity of predictions 
to values of nuisance parameters, and include correlations between rates and shapes, between signals and backgrounds, 
and between channels within experiments and between experiments. More discussion on this topic can be found in the 
individual analysis notes [ij through [ij]. Here we consider only the largest contributions and correlations between 
and within the two experiments. 

1. Correlated Systematics between CDF and DO 

The uncertainties on the measurements of the integrated luminosities are 6% (CDF) and 6.1% (DO). Of these 
values, 4% arises from the uncertainty on the inelastic pp scattering cross section, which is correlated between CDF 
and DO. CDF and DO also share the assumed values and uncertainties on the production cross sections for top- 
quark processes {tt and single top) and for electroweak processes {WW, WZ, and ZZ). In order to provide a 
consistent combination, the values of these cross sections assumed in each analysis are brought into agreement. We 
use CTff = 7.04j^Q'3g (scale) ± 0.14(PDF) ± 0.30(mass), following the calculation of Moch and Uwer [G^, assuming a 
top quark mass m* = 173.0 ± 1.2 GeV/c^ [24,], and using the MSTW2008nnlo PDF set [31]. Other calculations of att 
are similar [63 |. 

For single top, we use the NLL i-channel calculation of Kidonakis J6a], which has been updated using the 
MSTW2008nnlo PDF set [U [H]. For the s-channel process we use again based on the MSTW2008nnlo 

PDF set. Both of the cross section values below are the sum of the single t and single t cross sections, and both 
assume rrit = 173 ± 1.2 GeV. 

CTt-chan = 2.10 ± 0.027(scale) ± 0.18(PDF) ± 0.045(mass)pb. (5) 

CT^.^han = 1-046 ± 0.006(scale) ± 0.059 (PDF) ± 0.030 (mass) pb. (6) 

Other calculations of usingicTop are similar for our purposes [gsI] . 

mcfm [20] has been used to compute the NLO cross sections for WW, WZ, and ZZ production [69]. Using a scale 
choice Ho = + p'^{V) and the MSTW2008 PDF set the cross section for inclusive W+W~ production is 

aw+w- = ll-34tE!:l^ (scale) t;^;35(PDF)pb (7) 

and the cross section for inclusive W^Z production is 

aw±z = 3.22^0:20 (scale) tVol (PDF) pb (8) 

For the Z, leptonic decays are used in the definition of the cross section, and we assume both 7 and Z exchange. The 
W^Z cross section quoted above involves the requirement 75 < me+g~ < 105 GeV for the leptons from the neutral 
current exchange. The same dilepton invariant mass requirement is applied to both sets of leptons in determining the 
ZZ cross section which is 

azz = 1.20t°:l]l (scale) t°o'ot (PDF) pb (9) 

For the diboson cross section calculations, we use jry^] < 5 for all calculations. Loosening this requirement to include 
all leptons leads to ~-|-0.4% change in the predictions. Lowering the factorization and renormalization scales by a 
factor of two increases the cross section, and raising the scales by a factor of two decreases the cross section. The PDF 
uncertainty has the same fractional impact on the predicted cross section independent of the scale choice. All PDF 
uncertainties are computed as the quadrature sum of the twenty 68% C.L. eigenvectors provided with MSTW2008 
(MSTW2008nlo68cl). We furthermore constrain a^+w- in the signal regions of the channels in the process of setting 
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our limits, either by integration over the uncertain parameters, or by a direct fit, depending on the method. Our 
posterior constraint on (Jw+w- has a fractional uncertainty of ±4%, which includes the prior theoretical constraint 
of ±6%, indicating that the data and the theoretical prediction are approximately equally constraining on ctw+w-- 

In many analyses, the dominant background yields are calibrated with data control samples. Since the methods 
of measuring the multijet ("QCD") backgrounds differ between CDF and DO, and even between analyses within 
the collaborations, there is no correlation assumed between these rates. Similarly, the large uncertainties on the 
background rates for W+hea.vy flavor (HF) and Z+heavy flavor are considered at this time to be uncorrelated, as 
both CDF and DO estimate these rates using data control samples, but employ different techniques. The calibrations of 
fake leptons, unvetoed 7 — >■ e+e" conversions, 6-tag efficiencies and mistag rates are performed by each collaboration 
using independent data samples and methods, and are therefore also treated as uncorrelated. 

2. Correlated Systematic Uncertainties for CDF 

The dominant systematic uncertainties for the CDF analyses are shown in the Appendix in Tables IIX[ 1X1 and IXII 
for the H -> W+W' e^i'i"^ v cha nnels, in Table EI] for the WH WWW I'^l'^ and WH WWW 
^±£/±£//T channels, and in Table IXllIl for the ZH ZWW i^e^l'^ channels. Each source induces a correlated 
uncertainty across all CDF channels' signal and background contributions which are sensitive to that source. Shape 
dependencies of templates on jet energy scale and "FSR" ) are taken into account in the analyses (see tables) . For 
H — W^W~ , the largest uncertainties on signal acceptance originate from Monte Carlo modeling. Uncertainties on 
background event rates vary significantly for the different processes. The backgrounds with the largest systematic 
uncertainties are in general quite small. Such uncertainties are constrained by fits to the nuisance parameters, and 
they do not affect the result significantly. Because the largest background contributions are measured using data, 
these uncertainties are treated as uncorrelated for the iJ — > 66 channels. The differences in the resulting limits when 
treating the remaining uncertainties as either correlated or uncorrelated, is less than 5%. 

3. Correlated Systematic Uncertainties for DO 

The dominant systematic uncertainties for the DO analyses are shown in the Appendix, in Tables IXIVl IXVl IXVIl 
IXVIIi IXVIIIl and IXIXI Each source induces a correlated uncertainty across all DO channels sensitive to that source. 
Wherever appropriate the impact of systematic effects on both the rate and shape of the predicted signal and back- 
ground is included. For the H — >■ M^+VF^and VH — > li^t'^ + ^analyses, a significant source of uncertainty is the 
measured efficiencies for selecting leptons. Significant sources for all analyses are the uncertainties on the luminosity 
and the cross sections for the simulated backgrounds. For analyses involving jets the determination of the jet energy 
scale, jet resolution and the multijet background contribution are significant sources of uncertainty. All systematic 
uncertainties arising from the same source are taken to be correlated among the different backgrounds and between 
signal and background. 



VI. COMBINED RESULTS 

Before extracting the combined limits we study the distributions of the log-likelihood ratio (LLR) for different hy- 
potheses, to quantify the expected sensitivity across the mass range tested. Figure [5] displays the LLR distributions for 
the combined analyses as functions of vtlh- Included are the median of the LLR distributions for the background-only 
hypothesis (LLR;,), the signal-plus-background hypothesis (LLR^+f,), and the observed value for the data (LLRobs)- 
The shaded bands represent the one and two standard deviation (ct) departures for LLR;, centered on the median. 
Table |IV] lists the observed and expected LLR values shown in Figure El 
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FIG. 6: Distributions of the log-likelihood ratio (LLR) as a function of Higgs mass obtained with the CLs method for the 
combination of all CDF and DO analyses. 



These distributions can be interpreted as follows: The separation between the medians of the LLRb and LLR^+b 
distributions provides a measure of the discriminating power of the search. The sizes of the one- and two-cr LLR^ 
bands indicate the width of the LLRf, distribution, assuming no signal is truly present and only statistical fluctuations 
and systematic effects are present. The value of LLRobs relative to LLR^+b and LLR;, indicates whether the data 
distribution appears to resemble what we expect if a signal is present (i.e. closer to the LLR^+b distribution, which 
is negative by construction) or whether it resembles the background expectation more closely; the significance of any 
departures of LLRobs from LLR^ can be evaluated by the width of the LLR^ bands. 

Using the combination procedures outlined in Section III, we extract limits on SM Higgs boson production a x 
B{H — !• X) in pp collisions at y/s = 1.96 TeV for 130 < niH < 200 GeV/c^. To facilitate comparisons with the 
SM and to accommodate analyses with different degrees of sensitivity, we present our results in terms of the ratio of 
obtained limits to the SM Higgs boson production cross section, as a function of Higgs boson mass, for test masses 
for which both experiments have performed dedicated searches in different channels. A value of the combined limit 
ratio which is equal to or less than one indicates that that particular Higgs boson mass is excluded at the 95% C.L. 
A value less than one indicates that a Higgs boson of that mass is excluded with a smaller cross section than the SM 
prediction, and that the SM prediction is excluded with more certainty than 95% C.L. 

The combinations of results of each single experiment, as used in this Tevatron combination, yield the following 
ratios of 95% C.L. observed (expected) limits to the SM cross section: 0.92 (0.93) for CDF and 0.75 (0.92) for DO at 
niH — 165 GeV/c^. Both collaborations independently exclude ttih — 165 GeV/c^ at the 95% C.L. 

The ratios of the 95% C.L. expected and observed limit to the SM cross section are shown in Figure [7] for the 
combined CDF and DO analyses. The observed and median expected ratios are listed for the tested Higgs boson 
masses in Table IVl for niH < 150 GeV/c^, and in Table IVll for mn > 155 GeV/c^, as obtained by the Bayesian and 
the CLs methods. In the following summary we quote only the limits obtained with the Bayesian method, which 
was chosen a priori. It turns out that the Bayesian limits are slightly less stringent. The corresponding limits and 
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TABLE IV: Log-likelihood ratio (LLR) values for the combined CDF + DO Higgs boson search obtained using the CLs method. 
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expected limits obtained using the CLg method are shown alongside the Bayesian limits in the tables. We obtain 
the observed (expected) values of 1.31 (0.92) at tuh = 155 GeV/c^, 0.54 (0.65) at mu = 165 GeV/c^, 1.13 (0.85) at 
niH = 175 GeV/c2 and 1.49 (1.30) at mu = 185 GeV/c^. 



TABLE V: Ratios of median expected and observed 95% C.L. limit to the SM cross section for the combined CDF and DO 
analyses as a function of the Higgs boson mass in GeV/c'^, obtained with the Bayesian and with the CLs method. 
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TABLE VI: Ratios of median expected and observed 95% C.L. limit to the SM cross 
analyses as a function of the Higgs boson mass in GeV/c^, obtained with the Bayesian 


section for the combined CDF and DO 
and with the CLs method. 
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We show in Figure[8]and list in Tabic IVlIl the observed l-CLg and its expected distribution for the background-only 
hypothesis as a function of the Higgs boson mass. This is directly interpreted as the level of exclusion of our search 
using the CLg method. The region excluded at the 95% C.L. agrees very well with that obtained via the Bayesian 
calculation. 
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In addition, we provide in Figure [HK and listed in Table IVllH) the values for the observed l-CLg+b and its expected 
distribution as a function of mn- The value CLg+b is the p- value for the signal-plus-background hypothesis. These 
values can be used to obtain alternative upper limits that are more constraining compared to those obtained via the 
CLg method. In such a formulation, the power of the search is limited at a level chosen a priori to avoid setting limits 
when the background model grossly overpredicts the data or the data exhibit a large background-like fluctuation {e.g., 
limit at the -Icr background fluctuation level.). Within Figure |9l 95% C.L. power-constrained limits can be found 
at the points for which l-CLg+b exceeds 95%. The expected range of exclusion is ^40% larger using PCL than the 
Bayesian and CLs limits quoted here. We continue our convention of quoting Bayesian and CLg limits however. 

In summary, we combine CDF and DO results on SM Higgs boson searches, based on luminosities up to 8.2 fb^^. 
Compared to our previous combination, more data have been added to the existing channels, additional channels 
have been included, and analyses have been further optimized to gain sensitivity. We use the recommendation of the 
PDF4LHC working group for the central value of the parton distribution functions and uncertainties (s^l • We use the 
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FIG. 7: Observed and expected (median, for the background-only hypothesis) 95% C.L. upper limits on the ratios to the SM 
cross section, as functions of the Higgs boson mass for the combined CDF and DO analyses. The limits are expressed as a 
multiple of the SM prediction for test masses (every 5 GeV/c^) for which both experiments have performed dedicated searches 
in different channels. The points are joined by straight lines for better readability. The bands indicate the 68% and 95% 
probability regions where the limits can fluctuate, in the absence of signal. The limits displayed in this figure are obtained with 
the Bayesian calculation. 
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FIG. 8: The exclusion strength 1-CLs as a function of the Higgs boson mass (in steps of 5 GeV/c^), as obtained with CLs 
method for the combination of the CDF and DO analyses. 



highest-order calculations of gg — > H, WH, ZH, and VBF theoretical cross sections when comparing our limits to 
the SM predictions at high mass. We include consensus estimates of the theoretical uncertainties on these production 
cross sections and the decay branching fractions in the computations of our limits. 

The 95% C.L. upper limit on Higgs boson production is a factor of 0.54 times the SM cross section for a Higgs 
boson mass of mn =165 GeV/c^. Based on simulation, the corresponding median expected upper limit is 0.65 times 
the SM cross section. Standard Model branching ratios, calculated as functions of the Higgs boson mass, are assumed. 

We choose to use the intersections of piecewise linear interpolations of our observed and expected rate limits in order 
to quote ranges of Higgs boson masses that are excluded and that are expected to be excluded. The sensitivities of our 
searches to Higgs bosons are smooth functions of the Higgs boson mass and depend most strongly on the predicted 
cross sections and the decay branching ratios (the decay H — W^W~ is the dominant decay for the region of highest 
sensitivity) . The mass resolution of the channels is poor due to the presence of two highly energetic neutrinos in signal 
events. We therefore use the linear interpolations to extend the results from the 5 GeV/c^ mass grid investigated 
to points in between. This procedure yields higher expected and observed interpolated rate limits than if the full 
dependence of the cross section and branching ratio were included as well, since the latter produces limit curves that 
are concave upwards. We exclude in this way the region 158 < mu < 173 GeV/c^ at the the 95% C.L. The expected 
exclusion region, given the current sensitivity, is 153 < mn < 179 GeV/c^. The excluded region obtained by finding 
the intersections of the linear interpolations of the observed 1 — CLg curve shown in Figure [5] is slightly larger than 
that obtained with the Bayesian calculation. As previously stated, and following the procedure used in previous 
combinations [31 , we make the a priori choice to quote the exclusion region using the Bayesian calculation. 

The results presented in this paper significantly extend the individual limits of each collaboration and those obtained 
in our previous combination. The sensitivity of our combined search is sufficient to exclude a Higgs boson at high 
mass and is expected to grow substantially in the future as more data are added and further improvements are made 
to the analysis techniques. 
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TABLE VII: The observed and expected 1-CLs values as functions of mn, for the combined CDF and DO Higgs boson searches. 
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FIG. 9: The signal exclusion p- value l-CLg+b as a function of the Higgs boson mass (in steps of 5 GeV/c^) for the combination 
of the CDF and DO analyses. 
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TABLE VIII: The observed and expected l-CLg+b values as functions of mn, for the combined CDF and DO Higgs boson 

searches. 
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TABLE IX: Systematic uncertainties on the signal and background contributions for CDF's H — >■ W'^W~ — >■ i^i'^ channels 
with zero, one, and two or more associated jets. These channels are sensitive to gluon fusion production (all channels) and 
WH, ZH and VBF production. Systematic uncertainties axe listed by name (see the original references for a detailed explanation 
of their meaning and on how they axe derived) . Systematic uncertainties for H shown in this table are obtained for niH = 160 
GeV/c^. Uncertainties are relative, in percent, and are symmetric unless otherwise indicated. The uncertainties associated with 
the different background and signal processed arc correlated within individual jet categories unless otherwise noted. Boldface 
and italics indicate groups of uncertainties which are correlated with each other but not the others on the line. 
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CDF: H W+W- e^e^ with two or more associated jets channel relative uncertainties (%) 
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TABLE X: Systematic uncertainties on the signal and background contributions for CDF's low-Af^^ H — > W^W^ — > l^^l'^ 
channel with zero or one associated jets. This channel is sensitive to only gluon fusion production. Systematic uncertainties are 
listed by name (see the original references for a detailed explanation of their meaning and on how they are derived) . Systematic 
uncertainties for H shown in this table are obtained for niH = 160 GeV /(?. Uncertainties are relative, in percent, and are 
symmetric unless otherwise indicated. The uncertainties associated with the different background and signal processed are 
correlated within individual categories unless otherwise noted. In these special cases, the correlated uncertainties are shown in 
either italics or bold face text. 
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1.6 












2.1 


Higher-order Diagrams 




10.0 


10.0 


10.0 


10.0 






Jet Energy Scale 


1.0 


2.3 


2.0 


12.9 


6.4 


1.3 


2.4 


Conversion Modeling 












10.0 




Jet Fake Rates 












18.4 




Lepton ID Efficiencies 


3.8 


3.8 


3.8 


3.8 


3.8 




3.8 


Trigger Efficiencies 


2.0 


2.0 


2.0 


2.0 


2.0 




2.0 


Luminosity 


3.8 


3.8 


3.8 


3.8 


3.8 




3.8 


Luminosity Monitor 


4.4 


4.4 


4.4 


4.4 


4.4 




4.4 
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TABLE XI: Systematic uncertainties on the signal and background contributions for CDF's H — > W'^W~ — >■ e^r^ and 
H W'^W~ — > /x^r^ channels. These channels are sensitive to gluon fusion production, WH,ZH and VBF production. 
Systematic uncertainties are listed by name (see the original references for a detailed explanation of their meaning and on how 
they are derived). Systematic uncertainties for H shown in this table axe obtained for mu = 160 GeV/c^. Uncertainties are 
relative, in percent, and are symmetric unless otherwise indicated. The uncertainties associated with the different background 
and signal processed arc correlated within individual categories unless otherwise noted. In these special cases, the correlated 
uncertainties are shown in either italics or bold face text. 



CDF: H W+W- 



'T^ channel relative uncertainties ( ) 



Contribution 


WW 


WZ 


ZZ 


tt 


Z ^TT 


Z^ii 


W+jet 




gg^H 


WH 


ZH 


VBF 


Cross section 


6.0 


6.0 


6.0 


10.0 


5.0 


5.0 






10.3 


5 


5 


10 


Measured W cross-section 














12 












PDF Model 


1.6 


2.3 


3.2 


2.3 


2.7 


4.6 


2.2 


3.1 


2.5 


2.0 


1.9 


1.8 


Higher order diagrams 


10 


10 


10 


10 


10 


10 




10 




10 


10 


10 


Conversion modeling 
















10 










Trigger Efficiency 


0.5 


0.6 


0.6 


0.6 


0.7 


0.5 


0.6 


0.6 


0.5 


0.5 


0.6 


0.5 


Lepton ID Emciency 


0.4 


0.5 


0.5 


0.4 


0.4 


0.4 


0.5 


0.4 


0.4 


0.4 


0.4 


0.4 


T ID Efficiency 


1.0 


1.3 


1.9 


1.3 


2.1 






0.3 


2.8 


1.6 


1.7 


2.8 


Jet into T Fake rate 


5.8 


4.8 


2.0 


5.1 




0.1 


8.8 






4.2 


4.0 


0.4 


Lepton into r Fake rate 


0.2 


0.1 


0.6 


0.2 




2.3 




2.1 


0.15 


0.06 


0.15 


0.11 


W+jet scale 














1.6 












MC Run dependence 


2.6 


2.6 


2.6 








2.6 












Luminosity 


3.8 


3.8 


3.8 


3.8 


3.8 


3.8 


3.8 


3.8 


3.8 


3.8 


3.8 


3.8 


Luminosity Monitor 


4.4 


4.4 


4.4 


4.4 


4.4 


4.4 


4.4 


4.4 


4.4 


4.4 


4.4 


4.4 




CDF: 


H ^ 


W+W 


~ ^ n 


"^T^ channel relative uncertainties (%) 








Contribution 


WW 


WZ 


ZZ 


tt 


Z ^TT 


Z^U 


W+jet 


W-y 


gg^H 


WH 


ZH 


VBF 


Cross section 


6.0 


6.0 


6.0 


10.0 


5.0 


5.0 






10.4 


5 


5 


10 


Measured W cross-section 














12 












PDF Model 


1.5 


2.1 


2.9 


2.1 


2.5 


4.3 


2.0 


2.9 


2.6 


2.2 


2.0 


2.2 


Higher order diagrams 


10 


10 


10 


10 








11 




10 


10 


10 


Trigger Efficiency 


1.3 


0.7 


0.7 


1.1 


0.9 


1.3 


1.0 


1.0 


1.3 


1.3 


1.2 


1.3 


Lepton ID Efficiency 


1.1 


1.4 


1.4 


1.1 


1.2 


1.1 


1.4 


1.3 


1.0 


1.0 


1.0 


1.0 


T ID Efficiency 


1.0 


1.2 


1.4 


1.6 


1.9 








2.9 


1.6 


1.7 


2.8 


,Jct into T Fake rate 


5.8 


5.0 


4.4 


4.4 




0.2 


8.8 






4.5 


4.2 


0.4 


Lepton into r Fake rate 


0.06 


0.05 


0.09 


0.04 




1.9 




1.2 


0.04 


0.02 


0.02 


0.04 


W+jet scale 














1.4 












MC Run dependence 


3.0 


3.0 


3.0 








3.0 












Luminosity 


3.8 


3.8 


3.8 


3.8 


3.8 


3.8 


3.8 


3.8 


3.8 


3.8 


3.8 


3.8 


Luminosity Monitor 


4.4 


4.4 


4.4 


4.4 


4.4 


4.4 


4.4 


4.4 


4.4 


4.4 


4.4 


4.4 
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TABLE XII: Systematic uncertainties on the signal and background contributions for CDF's WH — >• WWW — >■ ^^^'^ channel 
with one or more associated jets and WH — > WWW — > I'^i'^l"^ channel. These channels are sensitive to only WH and 
ZH production. Systematic uncertainties are listed by name (see the original references for a detailed explanation of their 
meaning and on how they arc derived). Systematic uncertainties for H shown in this table are obtained for mn = 160 GeV/c^. 
Uncertainties are relative, in percent, and are symmetric unless otherwise indicated. The uncertainties associated with the 
different background and signal processed arc correlated within individual categories unless otherwise noted. In these special 
cases, the correlated uncertainties are shown in either italics or bold face text. 



CDF: WH WWW channel relative uncertainties (%) 



Contribution 


WW 


WZ 


ZZ 




tt 


DY W^ Vy-hjet 


WH 


ZH 


Cross Section 


6.0 


6.0 


6.0 




7.0 


5.0 


5.0 


5.0 


Scale (Acceptance) 


-6.1 
















PDF Model (Acceptance) 


5.1 
















Higher-order Diagrams 




10.0 


10.0 




10.0 


10.0 10.0 


10.0 


10.0 


Conversion Modeling 












10.0 






Jet Fake Rates 












38.5 






Jet Energy Scale 


-I4.O 


-3.9 


-2.8 




-0.6 


-77 -n.6 


-1.0 


-0.7 


Charge Mismeasurement Rate 


19.0 








19.0 


19.0 






Lepton ID Efficiencies 


3.8 


3.8 


3.8 




3.8 


3.8 


3.8 


3.8 


Trigger Efficiencies 


2.0 


2.0 


2.0 




2.0 


2.0 


2.0 


2.0 


Luminosity 


3.8 


3.8 


3.8 




3.8 


3.8 


3.8 


3.8 


Luminosity Monitor 


I'l 


I'l 










I'l 




CDF 


': WH 


WWW ^. 




^ channel relative uncertainties (%) 






Contribution 


WZ 


zz 




Zl 




tt Fakes 


WH 


ZH 


Cross Section 


6.0 


6.0 




10.0 




7.0 


5.0 


5.0 


Higher-order Diagrams 


10.0 


10.0 




15.0 




10.0 


10.0 


10.0 


Jet Energy Scale 








-2.7 










Jet Fake Rates 












25.6 






6- Jet Fake Rates 












27.3 






MC Run Dependence 








5.0 










Lepton ID Efficiencies 


5.0 


5.0 








5.0 


5.0 


5.0 


Trigger Efficiencies 


2.0 


2.0 








2.0 


2.0 


2.0 


Luminosity 


3.8 


3.8 








3.8 


3.8 


3.8 


Luminosity Monitor 


4.4 


4.4 








4.4 


4.4 


4.4 
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TABLE XIII: Systematic uncertainties on the signal and background contributions for CDF's ZH — >• ZWW — >■ i'^i^i''^ 
channels with 1 jet and 2 or more jets. These channels are sensitive to only WH and ZH production. Systematic uncertainties 
are listed by name (see the original references for a detailed explanation of their meaning and on how they are derived). 
Systematic uncertainties for H shown in this table axe obtained for niH = 160 GeV /<? . Uncertainties are relative, in percent, 
and are symmetric unless otherwise indicated. The uncertainties associated with the different background and signal processed 
arc correlated within individual categories unless otherwise noted. In these special cases, the correlated uncertainties are shown 
in either italics or bold face text. 



CDF: ZH ZWW t't^l'^ with one associated jet channel relative uncertainties (%) 



Contribution 


WZ 


ZZ 




tt 


Fakes WH 


ZH 


Cross Section 


6.0 


6.0 


10.0 


7.0 


5.0 


5.0 


Higher-order Diagrams 


10.0 


10.0 


15.0 


10.0 


10.0 


10.0 


Jet Energy Scale 


-7.6 


-2.3 


-5.3 


9.J, 


-9.0 


8.1 


Jet Fake Rates 










24.8 




6- Jet Fake Rates 








42.0 






MC Run Dependence 






5.0 








Lepton ID Efficiencies 


5.0 


5.0 




5.0 


5.0 


5.0 


Trigger Efficiencies 


2.0 


2.0 




2.0 


2.0 


2.0 


Luminosity 


3.8 


3.8 




3.8 


3.8 


3.8 


Luminosity Monitor 


44 


44 




44 


44 


4-4% 


CDF: ZH ZWW i^i^l'^ 


with two or 


more associated jets channel relative uncertainties (%) 




Contribution 


WZ 


ZZ 


Zl 


tt 


Fakes WH 


ZH 


Cross Section 


6.0 


6.0 


10.0 


7.0 


5.0 


5.0 


Higher-order Diagrams 


10.0 


10.0 


15.0 


10.0 


10.0 


10.0 


Jet Energy Scale 


-17.8 


-13.1 


-18.2 


-3.6 


-15.4 


-4.9 


Jet Fake Rates 










25.6 




&-Jet Fake Rates 








22.2 






MC Run Dependence 






5.0 








Lepton ID Efficiencies 


5.0 


5.0 




5.0 


5.0 


5.0 


Trigger Efficiencies 


2.0 


2.0 




2.0 


2.0 


2.0 


Luminosity 


3.8 


3.8 




3.8 


3.8 


3.8 


Luminosity Monitor 


44 


44 




44 


44 


44 
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TABLE XIV: Systematic uncertainties on the signal and background contributions for DO's H — >■ WW — ^ i^H!^ channels. 
Systematic uncertainties axe listed by name; see the original references for a detailed explanation of their meaning and on how 
they are derived. Shape uncertainties are labeled with the "s" designation. Systematic uncertainties shown in this table are 
obtained for the mn = 165 GeV/c^ Higgs selection. Uncertainties are relative, in percent, and are symmetric unless otherwise 
indicated. 

DO: H -J- WW f^i'^ channels relative uncertainties (%) 



Contribution 


Diboson 


U 


W +jet/'y 


tt 


Multijet 


H 


Luminosity/Normalization 


6 


6 


6 


6 


30 


6 


Cross Section 


7 


5 


6 


10 




7-33 


PDF 


2.5 


2.5 


2.5 


2.5 




8-30 


EM Identification 


2.5 


2.5 


2.5 


2.5 




2.5 


Muon Identification 


4 


4 


4 


4 




4 


Vertex Confirmation (s) 


2-6 


1-7 


1-6 


1-8 




1-8 


Jet identification (s) 


2-5 


2-5 


2-5 


2-5 




2-5 


Jet Energy Scale (s) 


2-3 


1-4 


1-8 


1-4 




1-10 


Jet Energy Resolution(s) 


1-4 


1-4 


1-12 


1-3 




1-12 


B-tagging 


10 


10 


10 


5 
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TABLE XV: Systematic uncertainties on the signal and background contributions for DO's H W~^W~ — ^ iivThadi' channel. 
Systematic uncertainties are listed by name; see the original references for a detailed explanation of their meaning and on 

how they arc derived. Shape uncertainties arc labeled with the shape designation (s). Systematic uncertainties shown in this 
table are obtained for the mn = 165 GeV/c^ Higgs selection. Uncertainties are relative, in percent, and are symmetric unless 
otherwise indicated. 

DO: H — ^ W^W' — ^ fii'Thadt^ channel relative uncertainties (%) 



Contribution 


Diboson 


Z/7* ^ U 


W + jets 


tt 


Multijet 


H 


Luminosity ((Tinei(pp)) 


4.6 


4.6 




4.6 




4.6 


Luminosity Monitor 


4.1 


4.1 




4.1 




4.1 


Trigger 


5.0 


5.0 




5.0 




5.0 


Lepton ID 


3.7 


3.7 




3.7 




3.7 


EM veto 


5.0 






5.0 




5.0 


Tau Energy Scale (s) 


1.0 


1.1 




<1 




<1 


Jet Energy Scale (s) 


8.0 


<1 




1.8 




2.5 


Jet identification (s) 


<1 


<1 




7.5 




5.0 


Multijet (s) 










20-50 




Cross Section 


7.0 


4.0 




10 




10 


Modeling 


1.0 




10 






3.0 
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TABLE XVI: Systematic uncertainties on the signal and background contributions for DO's WH — >• WWW i'^i'^ channel. 
Systematic uncertainties axe listed by name; see the original references for a detailed explanation of their meaning and on how 
they are derived. Shape uncertainties are labeled with the "shape" designation. Systematic uncertainties for signal shown 
in this table are obtained for tuh = 165 GeV/c^. Uncertainties are relative, in percent, and are symmetric unless otherwise 
indicated. 

DO: VH + X Run Ila channel relative uncertainties (%) 



Contribution 


WZ/ZZ 


W+jet 


ChargeFlip 


Multijet 


VH UX 


Cross section 


7 


6 











Normalization 


4 


4 











Trigger (/x/x) 














2 


LeptonID (ee) 


8.6 


8.6 








8.6 


LeptonID (///u) 


4 


4 








4 


LeptonID (e/n) 


6.3 


6.3 








6.3 


T j_TT\ / 

JetlD/JES 


2 


2 








2 


Jet-Lepton Fake 





20 











Instrumental (ee) 











52 


44 


Instrumental {efj, 














29 


Instrumental (/tt/tt) 











155 


42 


Instrumental Model 






shape 


shape 






DO: VH^i"''! 


+ X Run lib channel relative uncertainties (%) 




Contribution 


WZ/ZZ 


W+jet 


ChargeFlip 


Multijet 


VH eix 


Cross section 


7 


6 











Normalization 


4 


4 











Trigger (/x/x) 














5 


LeptonID (ee) 


8.6 


8.6 








8.6 


LeptonID (/i/i) 


4 


4 








4 


LeptonID (e/i) 


6.3 


6.3 








6.3 


JetlD/JES 


2 


2 








2 


Jet-Lepton Fake 





20 











Instrumental (ee) 











23 


31 


Instrumental (e/i) 














19 


Instrumental (nn) 











43 


28 


Instrumental Model 






shape 


shape 
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TABLE XVII: Systematic uncertainties on the signal and background contributions for DO's H — > WW* — > Ivjj electron and 
muon channels. Systematic uncertainties are listed by name; see the original references for a detailed explanation of their 
meaning and on how they are derived. Signal uncertainties are shown for tuh ~ 160 GeV/c^ for all channels except for WH, 
shown for niH = 115 GeV/c^. Those affecting the shape of the RF discriminant are indicated with "Y." Uncertainties are 
listed as relative changes in normalization, in percent, except for those also marked by "S," where the overall normalization is 
constant, and the value given denotes the maximum percentage change from nominal in any region of the distribution. 



DO: H — > WW* — > Ivjj Run II channel relative uncertainties (%) 



Contribution 


Shape 


W+jcts 


Z+jets 


Top 


Diboson 


99 


qq qqH 


WH 


Jet energy scale 


Y 




< 0.1 


±0.7 


±3.3 


,+5.7\ 
\-A.0) 


±1.5 




Jet identification 


Y 


±6.6^ 


< 0.1 


±0.5 


±3.8 


±1.0 


±1.1 


±1.0 


Jet resolution 


Y 


aif 


< 0.1 


±0.5 


1-0.5,' 


C+3.o^ 

1-0.5,' 


±0.8 


±1.0 


Association of jets with PV 


Y 


±3.2^ 


±1.3^ 


±1.2 


±3.2 


±2.9 


±2.4 


\-0.2) 


Luminosity 


N 


n/a 


n/a 


±6.1 


±6.1 


±6.1 


±6.1 


±6.1 


Muon trigger 


Y 


±0.4^ 


< 0.1 


< 0.1 


< 0.1 


< 0.1 


< 0.1 


< 0.1 


Electron identification 


N 


±4.0 


±4.0 


±4.0 


±4.0 


±4.0 


±4.0 


±4.0 


Muon identification 


N 


±4.0 


±4.0 


±4.0 


±4.0 


±4.0 


±4.0 


±4.0 


ALPGEN tuning 


Y 


±1.1^ 


±0.3^ 


n/a 


n/a 


n/a 


n/a 


n/a 


Cross Section 


N 


±6 


±6 


±10 


±7 


±10 


±10 


±6 


Heavy-fiavor fraction 


Y 


±20 


±20 


n/a 


n/a 


n/a 


n/a 


n/a 


PDF 


Y 


±2.0^ 


±0.7^ 


< 0.1^ 


< 0.1^ 


< O.l'^ 


< O.l'^ 


< 0.1^ 



Multijet Background 



Y 



Electron channel 
±6.5 



Muon channel 
±26 
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TABLE XVIII: Systematic uncertainties on the signal and background contributions for DO's rrjj Run lib charmel. Systematic 
uncertainties for the Higgs signal shown in this table are obtained for niH = 135 GeV/c^. Systematic uncertainties are listed 
by name; see the original references for a detailed explanation of their meaning and on how they are derived. Uncertainties are 
relative, in percent, and are symmetric unless otherwise indicated. A systematic is denoted as flat if it affects the normalization 
only, and as 'shape' otherwise. 

DO: HThadjj Run lib channel relative uncertainties (%) 



Contribution VH Signal VBF Signal GGF Signal W + jets Z + jets Top diboson Multijet 



Luminosity (DO specific) 


4.1 


4.1 


4.1 


4.1 


4.1 


4.1 


4.1 




Luminosity (Tevatron common) 


4.6 


4.6 


4.6 


4.6 


4.6 


4.6 


4.6 




At ID 


2.9 


2.9 


2.9 


2.9 


2.9 


2.9 


2.9 




H trigger 


8.6 


8.6 


8.6 


8.6 


8.6 


8.6 


8.6 




r energy correction 


9.8 


9.8 


9.8 


9.8 


9.8 


9.8 


9.8 




r track efficiency 


1.4 


1.4 


1.4 


1.4 


1.4 


1.4 


1.4 




r selection by type 


12,4.2,7 


12,4.2,7 


12,4.2,7 


12,4.2,7 


12,4.2,7 


12,4.2,7 


12,4.2,7 




Cross section 


6.2 


4.9 


33 


6.0 


6.0 


10.0 


7.0 




GGF Signal PDF 






29 












GGF HpT Reweighting (Shape) 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 




Vertex confirmation for jets 


4.0 


4.0 


4.0 


4.0 


4.0 


4.0 


4.0 




Jet ID (Shape) 


~10 


~10 


~10 


~10 


~10 


~10 


~10 




Jet Energy Resolution (Shape) 


~10 


~10 


~10 


~10 


-10 


~10 


~10 




Jet energy Scale (Shape) 


~15 


~15 


~15 


~15 


~15 


~15 


~15 




Jet pT 


5.5 


5.5 


5.5 


5.5 


5.5 


5.5 


5.5 




PDF reweighting 


2 


2 


2 


2 


2 


2 


2 





Multijet Normalization - - - .... 5.3 

Multijet Shape - - - .... 



DO: eThadjj R-un lib relative uncertainties (%) 



Contribution VH Signal VBF Signal GGF Signal W + jets Z + jets Top diboson Multijet 



Luminosity (DO specific) 


4.1 


4.1 


4.1 


4.1 


4.1 


4.1 


4.1 


Luminosity (Tevatron common) 


4.6 


4.6 


4.6 


4.6 


4.6 


4.6 


4.6 


EM ID 


4 


4 


4 


4 


4 


4 


4 


e trigger 


2 


2 


2 


2 


2 


2 


2 


r energy correction 


9.8 


9.8 


9.8 


9.8 


9.8 


9.8 


9.8 


r track efficiency 


1.4 


1.4 


1.4 


1.4 


1.4 


1.4 


1.4 


r selection by type 


12,4.2,7 


12,4.2,7 


12,4.2,7 


12,4.2,7 


12,4.2,7 


12,4.2,7 


12,4.2,7 


Cross section 


6.2 


4.9 


33 


6.0 


6.0 


10.0 


7.0 


GGF Signal PDF 






29 










GGF HpT Reweighting (Shape) 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


Vertex confirmation for jets 


4.0 


4.0 


4.0 


4.0 


4.0 


4.0 


4.0 


Jet ID (Shape) 


~10 


~10 


~10 


~10 


~10 


~10 


~10 


Jet Energy Resolution (Shape) 


~10 


~10 


~10 


~10 


~10 


~10 


~10 


Jet energy Scale (Shape) 


~15 


~15 


~15 


~15 


~15 


~15 


~15 


Jet pT 


5.5 


5.5 


5.5 


5.5 


5.5 


5.5 


5.5 


PDF reweighting 


2 


2 


2 


2 


2 


2 


2 



Multijet Normalization - - - .... 4.7 

Multijet Shape - - - .... 
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TABLE XIX: Systematic uncertainties on the signal and background contributions for DO's H ^ channel. Systematic 
uncertainties for the Higgs signal shown in this table are obtained for thh = 125 GeV/c^. Systematic uncertainties are listed 
by name; see the original references for a detailed explanation of their meaning and on how they are derived. Uncertainties are 
relative, in percent, and are symmetric unless otherwise indicated. 

DO: — >■ 77 channel relative uncertainties (%) 



Contribution Background Signal 

Luminosity 6 6 

Acceptance - 2 

electron ID efficiency 2 - 

electron track-match inefficiency 10 - 

Photon ID efficiency 3 3 

Photon energy scale 2 1 

Cross Section 4 10 

Background subtraction 15 



